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ABSTRACT 

The O I 135.56nm line is covered by NASA’s Interface Region Imaging Spectrograph (IRIS) small explorer 
mission which studies how the solar atmosphere is energized. We here study the formation and diagnostic 
potential of this line by means of non-LTE modelling employing both ID semi-empirical and 3D radiation- 
Magneto Hydrodynamic (RMHD) models. We study the basic formation mechanisms and derive a quintessen¬ 
tial model atom that incorporates the essential atomic physics for the formation of the 01135.56nm line. This 
atomic model has 16 levels and describes recombination cascades through highly excited levels by effective 
recombination rates. The ionization balance O I/O II is set by the hydrogen ionization balance through charge 
exchange reactions. The emission in the O I 135.56nm line is dominated by a recombination cascade and the 
line is optically thin. The Doppler shift of the maximum emission correlates strongly with the vertical velocity 
in its line forming region, which is typically located at 1.0 - 1.5 Mm height. The total intensity of the line emis¬ 
sion is correlated with the square of the electron density. Since the 01135.56 nm line is optically thin, the width 
of the emission line is a very good diagnostic of non-thermal velocities. We conclude that the 01135.56 nm line 
is an excellent probe of the middle chromosphere, and compliments other powerful chromospheric diagnostics 
of IRIS such as the Mg II h & k lines and the C II lines around 133.5 nm. 

Subject headings: line: formation - radiative transfer - Sun: atmosphere - Sun: chromosphere 


1. INTRODUCTION 

The Ol 135.56 nm intersystem line (2s^2p^ 3s^S2— 
2s^ 2p^ ^P 2 ) is covered by the NASA/SMEX mission Interface 
Region Imaging Spectrograph (IRIS, De Pontieu et al. 2014). 
Unlike other strong lines in the IRIS spectral passbands, such 
as the h & k resonance lines of singly ionised magnesium 
(Leenaarts et al. 2013a,b; Pereira et al. 2013) or the C II 
133.5 nm multiplet (Rathore & Carlsson 2015; Rathore et al. 
2015) that form in the upper chromosphere to lower transi¬ 
tion region, the O I 135.56 nm line generally forms at mid- 
chromospheric heights, see below. Hence developing its po¬ 
tential diagnostics will provide a useful tool to exploit the data 
from IRIS and map different regions in the atmosphere. 

Already the SKYLAB mission discovered some interesting 
behaviour of the OI 135.56 nm line and its neighbouring per¬ 
mitted transition of CI during a solar flare (Cheng et al. 1980). 
The C I/O I line ratio got greatly enhanced during the solar 
flare compared with its value in the quiet sun. The lack of 
a theoretical explanation makes it hard to interpret what the 
mechanism is behind such a dramatic change. Cheng et al. 
(1980) suggested that this line ratio can be sensitive to the 
chromospheric electron density, a critical parameter for un¬ 
derstanding the energy balance of chromospheric plasmas. 

Non-Thermodynamic equilibrium (non-LTE) modelling of 
the neutral oxygen allowed resonance lines at 130.22, 130.49 
and 130.60 nm in the triplet system and the intercombination 
lines at 135.56 and 135.85 nm was carried out by Skelton & 
Shine (1982). The hydrogen Lyj3 line is very close in wave¬ 
length to an oxygen resonance line in the triplet system and ra¬ 
diative excitation of the oxygen line by hydrogen Lyj3 photons 
was found to be an important excitation mechanism for the OI 
130nm resonance lines. The excitation by Lyj3 pumping is 
transferred to the lower excited states in the triplet system by 
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radiative cascading. Skelton & Shine (1982) also confirmed 
that the ionization of neutral oxygen is dominated by charge- 
exchange reactions with hydrogen and therefore is determined 
by the ionization degree of hydrogen. They only included one 
state (2s^ 2p^ 3s ^S 2 ) in the quintet system, which is the upper 
level of the intercombination lines. Without the higher excited 
states in the quintet system, the authors concluded that the up¬ 
per level of the O I 135.56 nm line was mainly populated by 
collisional excitation. They also reported that the formation of 
the O I intersystem lines was less optically thick than the O I 
resonance lines at 130 nm, due to the lack of a central reversal 
in the line profile. 

The basic formation of the neutral oxygen lines was revis¬ 
ited by Carlsson & Judge (1993) with applications to a range 
of cool stars and solar chromospheric models, with a particu¬ 
lar focus on the resonance lines at 130nm. The authors con¬ 
firmed that the ionization degree of neutral oxygen is set by 
charge transfer with hydrogen, and the formation of the res¬ 
onance lines is dominated by the Lyj8 pumping effect. The 
model atom developed by Carlsson & Judge (1993) was de¬ 
veloped for the purpose of studying the resonance lines in the 
triplet system. Since the ionization balance is determined by 
charge transfer reactions with hydrogen and the excitation is 
dominated by hydrogen Lyj8 pumping, highly excited states 
do not have to be included when the interest is focussed on the 
resonance triplet. Their model atom has more excited states 
in the quintet system compared to Skelton & Shine (1982) but 
the model atom is not complete enough for the study of the 
intersystem lines that may be populated through a recombina¬ 
tion cascade in the quintet system. 

Fabbian et al. (2009) studied a number of O I lines for ele¬ 
mental abundance analysis, among them the 01777 nm lines, 
which belong to the quintet system. The upper levels of the 
777 nm lines are populated mostly through radiative cascading 
processes within the quintet system, therefore such high-lying 
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excited states were included in their study of the non-LTE ef¬ 
fects on the oxygen abundance. There was, however, no Lyj8 
pumping effects included in their modelling. Due to their in¬ 
clusion of the highly excited levels in the quintet system, we 
use their model atom as a basis for constructing our simplified 
model atom. 

The layout of this paper is as follows: In Section 2 we 
describe our radiative transfer computations and model at¬ 
mospheres. In Section 3 we construct a simplified model 
atom for the study of the O I 135.56nm line and describe this 
quintessential model atom. In Sections 4.1 and 4.2 we discuss 
the basic formation mechanism and the effects from the hy¬ 
drogen solution based on a ID semi-empirical atmosphere. In 
Section 4.3 we discuss different formation scenarios in a 3D 
atmosphere, in particular how velocity fields affect the line 
profile.. In Section 5 we present syntethic spectra and the po¬ 
tential diagnostics of the O I 135.56nm line. In Section 6 we 
show the comparison of our simulation with observations, and 
we conclude in Section 7. 

2. METHOD 

2.1. Radiative transfer computations 

We use RH (Uitenbroek 2001) to solve the non-LTE radia¬ 
tive transfer problem. RH is a multilevel accelerated lambda 
iteration code for radiative transfer calculations that includes 
partial frequency redistribution (PRD) treatment of the line 
profiles and also blending between lines. We use the RH ID 
version to study the basic formation mechanism in Section 4 
with the EALC (Eontenla et al. 1993) atmosphere. Eor syn¬ 
thetic profiles from the 3D atmosphere, we solve the radiative 
transfer column-by-column as a ID problem using the L5D 
version of RH (Pereira & Uitenbroek 2015). We take the non- 
LTE population densities of hydrogen from the model atmo¬ 
spheres and solve for the non-LTE population densities of sil¬ 
icon (the element responsible for the dominant background 
opacity) in addition to oxygen. 

2.2. Model atmospheres 

We use two model atmospheres for our study. In Section 4 
we use the semi-empirical solar chromosphere model EALC 
(Eontenla et al. 1993) to study the basic formation mecha¬ 
nisms. The EALC model is a one-dimensional atmosphere in 
hydrostatic equilibrium that describes the atmosphere in an 
ill-determined averaged fashion. It is thus far from catch¬ 
ing the dynamic state of the solar chromosphere but serves 
the purpose of providing a test-ground with reasonable atmo¬ 
spheric parameters for studying atomic processes to develop 
an understanding of the basic formation mechanisms. This 
understanding is essential for the development of a simpli¬ 
fied atomic model that encompasses the essentials while be¬ 
ing small enough to permit radiative transfer solutions in a 3D 
atmosphere. 

To study the formation properties in a more realistic at¬ 
mospheric model that catches some of the dynamic state 
of the solar chromosphere, we use a snapshot calculated 
with the 3D Radiation Magnetohydrodynamic (RMHD) code 
Bifrost (Gudiksen et al. 2011). We use a snapshot that ex¬ 
tends 24x24x16.8 Mm^ in physical space, with resolution of 
504x504x496 grid-points. It covers the upper convection zone 
(from 2.5 Mm below the height where optical depth is unity 
at 500 nm, the zero point of our height scale), the chromo¬ 
sphere, transition region and lower corona and includes mag¬ 
netic field with two main patches of opposite polarity sepa¬ 
rated by 8 Mm. The average unsigned field strength is 50G. 


The simulation cube is snapshot 385 of the en024048Jiion 
simulation that is available from the European Hinode Sci¬ 
ence Data Centre (http://www.sdc.uio.no/searcti/simulations). 
A detailed description of the simulation en024048Jiion is 
given in Carlsson et al. (2015). This is the same simulation 
cube that has been used in the previous papers on the forma¬ 
tion of IRIS diagnostics (Leenaarts et al. 2013a,b; Pereira et al. 
2013, 2015) and a number of other papers on line formation 
under solar chromospheric conditions (Leenaarts et al. 2012; 
Stepan et al. 2012; de la Cruz Rodriguez et al. 2013). 

3. QUINTESSENTIAL MODEL ATOM 

As a starting point for developing a simplified atomic model 
we use the 54-level model atom of Eabbian et al. (2009). Eor 
the source of the atomic data, see that reference. This 54-level 
model atom contains extra high lying excited states and up¬ 
dated collisional data compared with the 14-level model atom 
in Carlsson & Judge (1993). We simplify the 54 level atom to 
a model atom with 16 levels to meet our purpose. We do this 
by replacing all levels above the most highly excited level in 
Carlsson & Judge (1993) by effective recombination rates that 
catches the recombination cascade through these levels with¬ 
out solving for their population density, see below. Our re¬ 
sulting model atom is with 16 levels instead of 14 levels as the 
2p^3p ^Pi,2p^3p ^P 2 and 2p^3p^P3 are treated as separated 
levels, instead of a single merged term as in Carlsson & Judge 
(1993). A major uncertainty in the model atom is the impor¬ 
tance of collisions with neutral hydrogen. These collisions 
may couple the quintet system with the triplet system and 
make the O I 135.56nm line emission sensitive to the pump¬ 
ing by hydrogen Lyj3 radiation; a process that takes place in 
the triplet system. Collisional data for collisions with neu¬ 
tral hydrogen are largely unknown. A common approach in 
stellar abundance studies is to use a simplified recipe (Drawin 
1968, 1969) together with a scaling factor determined from 
a fit to observations. This approach is not ideal since a mis¬ 
match between observations and the non-LTE solution may 
come from the deficiency of other approximations (e.g., ID 
plane-parallel, hydrostatic equilibrium atmospheres) and not 
from lack of neutral hydrogen collisions. Eabbian et al. (2009) 
included the Drawin recipe in their calculation but did not in¬ 
clude the effect of Lyj3 pumping. The effect on the 01777 nm 
lines was rather small. We have made experiments using the 
ID EALC model and find that neutral hydrogen collisions will 
make the Lyj3 pumping contribute to the population of the up¬ 
per level of the O I 135.56 nm line and increase the intensity 
by up to 50%. However, quantum mechanical calculations of 
neutral hydrogen collisional cross-sections for other elements 
(none have been made for O I yet) indicate that the Drawin 
recipe may overestimate the rates by several orders of magni¬ 
tude (see Caccin et al. 1993, for a discussion). Eor this reason 
we have chosen not to include collisions with neutral hydro¬ 
gen but in light of the dependency we have found, it is im¬ 
portant to recognize this source of uncertainty in the resulting 
intensities. 

Highly excited levels are important for providing recombi¬ 
nation channels from the continuum. After the recombina¬ 
tion, we expect further downward cascading through allowed 
radiative transitions. In the chromosphere, we expect all these 
transitions to be optically thin such that we can neglect the 
reverse transitions (photoionization and radiative excitation). 
This makes it possible to introduce effective recombination 
rates that take into account the recombination channels with¬ 
out treating the levels in detail. 
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In general, we have a radiative recombination coefficient 
given by 


Table 1 

Atomic parameters for O. 


Rci = ^nNe^{ 
2gc 


InniekT 



Oic{v) , 2hv^ 

hv 


)e kTdv 


( 1 ) 

where the index c stands for the continuum, the index i the 
level it recombines to, Ne the electron density, gc, gi the sta¬ 
tistical weights. Me the electron mass, T the temperature, h 
the Planck constant, k the Boltzmann constant, Xd the ioniza¬ 
tion energy, v frequency, Vq the threshold frequency, Gic the 
photoionization crossection, and c the speed of light. 

If we factor out the electron density, Ne, the rest of the ex¬ 
pression in Eq. (1) will be just a function of temperature after 
carrying out the integration over frequency. We can thus ex¬ 
press the recombination to level i as: 


Rci=NeKi{T) (2) 

The advantage of such an approach is that if the net rate be¬ 
tween the continuum and the state i can be well approximated 
as this recombination process (meaning that the reverse pho¬ 
toionization is negligible), then we can replace the radiative 
recombination with a fixed rate that is only a function of tem¬ 
perature. If the electrons in state i further cascade down to 
other states with allowed transitions, e.g., to states j and k, 
then the recombination to state i can be split into recombina¬ 
tion processes between the continuum and the states j and k, 
with a split given by the transition probabilities for sponta¬ 
neous deexcitation, Aij and 


Ki{T) = T^hrKiiT) + - 


Mj 


Mi -^Mk 


KiiJ) 


(3) 




where • is the contribution to from recombination 
through level i. This process can be repeated for all the levels 
where we want to take into account the recombination cas¬ 
cade without including the level in the model atom. For each 
level we keep in the model atom, we thus end up with recom¬ 
bination contributions from all levels we exclude that couple 
to our level through allowed transitions. These contributions 
are just summed up to get the total effective recombination to 
that level. 

Figure 1 shows the term diagram of our quintessential 
model atom with 16 levels. The level energies and designa¬ 
tions are given in Table 1 . The effective recombination coeffi¬ 
cients to all levels are tabulated in Table 2 and we show them 
in Figure 2. 


4. BASIC FORMATION MECHANISM 
4.1. Ionization balance and the population of the upper level 

The ionization degree of the neutral oxygen is mainly de¬ 
cided collisionally through charge transfer with H and H+ 
(Judge 1986), 


0 + H+t^0++H (4) 

and the relation between the oxygen populations and the 
hydrogen populations follows Eq. (3) in Carlsson & Judge 
(1993): 

^(01) ^9 ^(HI) 

^(OII) 8^(HII)' ^ ^ 


Level 

Energy [cm 

Designation 

0 

0.000 

012 / ¥2 

1 

158.265 

012/^1 

2 

226.977 

Ol 2 /¥o 

3 

15867.862 

0 1 2/ ‘Z)2 

4 

33792.583 

012/ ‘So 

5 

73768.200 

012 / 733^550 

6 

76794.978 

Ol2p^3s^S\ 

7 

86625.757 

Ollp^ip^Pl 

8 

86627.778 

Ol2/j33/j5p2 

9 

86631.454 

Ol2/j33/j5p3 

10 

88630.977 

0\2p^2p^P 

11 

95476.728 

Ol2p^4s^Sl 

12 

96225.049 

Ol2p^4s^S\ 

13 

97420.748 

Ol2p^2d^ir 

14 

97488.476 

Ol2p^3d^D‘‘ 

15 

109837.020 

0 II ground 


The excited states in our atomic model belong to two main 
systems, the triplet and the quintet system. The upper level of 
the 135.56 nm line, 2p^ 3s ^82, belongs to the quintet system. 
The lower level of this line, 2p'^^P2, which is in the ground 
term, belongs to the triplet system. The 135.56 nm line is 
hence an intersystem line with low oscillation strength and 
therefore low opacity. The population and depopulation of 
the upper level is shown in Figure 3. The dominant chan¬ 
nel for electrons to come into the 2p^ 3s^S2 level is from the 
three levels in the 2p^3p^P term, (levels 7,8,9). The sec¬ 
ondary channel of incoming electrons is through collisions 
from the 2p^ 3s^S^ term (level 6), which is part of the triplet 
system. The 2p^3s^S2 level is depopulated by the emission 
lines at 135.56nm and 135.85nm, which correspond to radia¬ 
tive rates R5-0, R5-1. There is also a collisional rate, C5-0, 
corresponding to the 135.56 nm line, but its effect is rather 
small compared to the other channels. 

The upper level of the O I 135.56 nm line is thus populated 
mostly by radiative cascading from the continuum through the 
2p^3p^P term and is depopulated through the 135.56nm and 
the 135.85 nm line. As we will see in Section 5.1, due to the 
low oscillator strength of our line, it is usually optically thin. 

4.2. Effects of hydrogen 

The effects from hydrogen is twofold: one from the ioniza¬ 
tion degree of hydrogen, and the other from the radiation field 
of the Fyj3 line. 

As a primary channel, the states in the quintet system will 
be populated by a cascading process from the continuum and 
their population will therefore depend on the oxygen ioniza¬ 
tion degree. Since the oxygen ionization is dominated by 
charge transfer with neutral hydrogen and protons, higher hy¬ 
drogen ionization degree will result in higher oxygen ioniza¬ 
tion degree (see Equation 5). 

As a secondary channel the quintet levels can also be popu¬ 
lated through collisions from the triplet levels (C6-5 in Fig. 3), 
and with such a channel the Fyj8 radiation might have an 
effect on the 135.56nm emission as well. The reason is 
that the Fyj3 pumping of the O I 102.5 nm resonance line 
dominates the excitations in the triplet system (Carlsson & 
Judge 1993). We have made tests with and without the ef¬ 
fect of Fyj8 pumping included and find negligible effect on 
the O I 135.56 nm line if only collisions with electrons cou- 
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Figure 1. Term diagram of the 16 level OI quintessential model atom used in this study. The energy and designation of each numbered level is listed in table 1. 
The arrows show the direction of the net rates for each transition in the chromosphere. Vacuum wavelengths are shown in nm. The energies of the levels of the 
terms Ip^^P (levels 0,1,2) and 2p^ 3p^P (levels 7,8,9) are very close to each other and the energy differences within these terms are exaggerated in this diagram. 
The transitions to or from levels 7,8 and 9 are represented by the transition to and from level 8 {2p^ 3p^Pi) only. The O 1 135.56nm line that is the topic of this 
paper is shown in thick red, the 102.576 nm line that is pumped by the hydrogen Lyj8 line is in thick blue. 


Table 2 

Effective Recombination Coefficients 


Ki 

4500 

5160 

6370 

Temperature [K] 
7970 9983 

20420 

41180 

60170 

100000 

[10-28 ^ 3 ^- 1 ] 

254.126 

227.679 

190.951 

156.831 

126.345 

57.9456 

23.8875 

14.3800 

7.03006 

Ki3 [10-22 

45.3115 

40.6429 

34.0479 

27.9286 

22.5443 

10.3130 

4.26601 

2.55122 

1.25207 

R*12 [10-29 

164.436 

145.251 

118.711 

94.9460 

74.9470 

32.6631 

13.1234 

7.80194 

3.80117 

R*11 [10-28 

247.508 

221.465 

185.268 

151.301 

121.811 

55.5946 

22.9252 

13.7275 

6.72938 

Kw [10“^® m^s~^] 

85.4289 

76.4623 

63.7694 

52.0674 

41.9316 

19.1376 

7.90372 

4.72786 

2.32700 

Rig [10-28 

82.8493 

74.1332 

62.0482 

50.8526 

41.0086 

18.8720 

7.82051 

4.69823 

2.30519 

R*8 [10-28 

54.3390 

48.6406 

40.7773 

33.3889 

26.8978 

12.3791 

5.13489 

3.08086 

1.51400 

R*7 [10-29 m2i-‘] 

289.806 

259.521 

217.337 

177.804 

143.389 

65.9316 

27.3822 

16.4199 

8.07420 

[10-20 m2i-‘] 

376.007 

332.290 

272.020 

217.439 

171.643 

74.8017 

30.0632 

17.8556 

8.70174 

[10-29 

253.012 

226.401 

189.665 

154.888 

124.757 

57.0020 

23.5347 

14.0899 

6.90617 

R*2 [10-29 m2i-‘] 

260.786 

233.351 

195.260 

159.762 

128.709 

59.2080 

24.5266 

14.7083 

7.23052 

[10-28 

78.3535 

70.0987 

58.6938 

47.9973 

38.7065 

17.7829 

7.36847 

4.42349 

2.17186 

Rln [10-28 

130.846 

117.275 

98.0652 

80.2916 

64.5800 

29.7279 

12.3206 

7.39798 

3.62685 


pie the triplet and quintet system. Collisions with neutral hy¬ 
drogen may increase the coupling, but the collisional cross- 
sections are unknown and likely not large enough to affect the 
O I 135.56nm line, see Section 3 for a discussion. The main 
ionization/excitation and de-excitation channels are summa¬ 
rized in Figure 4. 

4.3. Line Formation in a 3D atmosphere 

We here illustrate the formation of the 01135.56 nm line by 
showing the formation in detail with the help of the four-panel 
diagrams introduced by Carlsson & Stein (1997). 


The emergent intensity along the normal in a ID plane par¬ 
allel semi-infinite atmosphere can be written as 

ly = f Sye ( 6 ) 

Jo 

where the source function (^v), opacity (Xv) and optical 
depth (Tv) are functions of frequency (v) and geometrical 
height (z). The integrand describes the local creation of pho¬ 
tons (SyXydz) and the fraction of those that escape (e~^^). We 
therefore define the contribution function to intensity on a ge- 
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Temperature [K] 

Figure 2. Effective recombination coefficients, see Table 2. The labels of 
the excited states are on top of each line with the same color. If two lines are 
very close, the label of the line with slightly lower value is shifted toward the 
right. 



Figure 3. Net rates to 2p^ 3s ^82, (level 5 in the 01 model atom (Figure 1)). 
Positive values denote a net rate coming into the level, negative values a rate 
going out of the level. Letters R and C in the label stands for radiative net 
rate and collisional net rate, respectively. The collisional net rates are also 
overlaid with star symbols. The dominant channel into the 2p^ 3s^82 level is 
from the three levels in the 2p^3p^P term, (levels 7,8,9, solid black). The 
secondary channel is through collisions from the 2p^3s^8^^ term (level 6, 
green dot-dot-dashed with stars). The 2p^ 3s^82 level is depopulated by the 
emission lines at 135.56nm (solid red) and 135.85nm (dashed yellow). There 
is also a collisional rate, C5-0, corresponding to the 135.56 nm line (dot-dash 
blue with stars). 



Figure 4. O I term diagram that summarizes the basic formation of the O I 
135.56 nm line (Section 4) . The solid lines represent radiative cascading, 
the dashed line radiative excitation and the dotted line stands for collisions. 
The ionization degree of O I is set by the ionization degree of hydrogen 
through charge transfer. The main rate into the quintet system, where the 
state 2p^2s^82 belongs, is through radiative cascading. Within the triplet 
system, the Lyj8 pumping drives photo-excitation which is followed by ra¬ 
diative cascading within the triplet system. There may, however, also be a 
channel from the triplet system into the quintet system through collisions, 
thus making the 01 135.56nm line sensitive to Lyj8 pumping. 

ometrical height scale as 

C,,{z)=Sye-^^Xy (7) 

Following Carlsson & Stein (1997) we can rewrite the con¬ 
tribution function as 

C/,{z)=5vTve-^''^. (8) 

Tv 

where the term has a maximum at Ty = 1 and repre¬ 

sents the Eddington-Barbier part of the contribution function, 
gives the source function contribution and the final term, 
^ picks out effects of velocity gradients in the atmosphere. 

Figure 5 shows the formation at an atmospheric column 
with moderate velocities. The continuum is formed at T = 1 
located at a height of 0.73 Mm where the temperature is 3 kK 
but with a decoupled source function at a radiation tempera¬ 
ture of 4 kK. The line core has optical depth unity only slightly 
higher, at z = 0.75 Mm, but the intensity is formed much 
higher up, at z = 1.2 Mm (shown by the contribution func¬ 
tion distribution in the lower right panel of Figure 5). The 
intensity comes from there because of the enhanced source 
function at this height (upper right panel). The high source 
function comes from a peak in the electron density (upper left 
panel). The line formation is thus optically thin in this case. 
The profile is not much modified by the velocity field and 
has a Gaussian shape with an 1/e width of 2.67 km s“^ which 
should be compared with the thermal width of 2.55 km s“^ 
for the atmospheric temperature of 6.3 kK at the formation 
height. 

Figure 6 shows the formation at a location where there is 
a substantial velocity field in the line forming region. The 
continuum is formed at a height of 0.7 Mm where the radia¬ 
tion temperature of the source function is 4.5 kK. The contin¬ 
uum intensity is therefore higher than in the case illustrated in 
Figure 5. The contribution function to intensity is significant 
over a large part of the atmosphere with a low optical depth. 
The line is thus optically thin. The line profile gets broad- 
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Figure 5. Intensity formation for the O I 135.56nm line in a column of the 
3D atmosphere characterised by moderate velocities. The label in the top-left 
corner of each panel gives the quantity shown as a grey-scale image as func¬ 
tion of wavelength (given as Doppler shift from the rest frequency, positive 
means redshift) and height z. The T = 1 height is given as a red dashed line 
and the vertical velocity as a blue solid line in all four panels with a white 
straight line denoting the rest frequency in the upper left panel. The green 
solid line in the upper left panel gives the electron density with a scale at the 
top. In the upper right panel we show the temperature as a yellow dashed line 
and the source function at a frequency given by the velocity profile as a green 
line, both with the scale given at the top. In the lower right panel we show 
the full contribution function (the product of the terms shown as grey-scale 
images in the other three panels) together with the emergent intensity profile 
with a scale to the right. 


Ne[10'"m"^] T[kK] 

1.6 1.4 1.2 1.0 0.8 0.6 0.4 6.5 6.0 5.5 5.0 4.5 4.0 3.5 
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Figure 6. As Figure 5 but for a column with a large velocity amplitude 
through the atmosphere(~ 15 km s“^). The line profile is highly asymmetric 
and the line is broadened by the velocity range within the O I line forming 
region (0.7-1.4 Mm in this case). Note the expanded wavelength scale com¬ 
pared with Figure 5. 


ened by the velocity gradients in the line forming region, in 
the 0.7-1.4 Mm height range. The 1/e width of the profile 
is 9.8 km s“^ substantially larger than the thermal width at 
the temperature of the line forming region. This line-width 
- velocity profile relation can be used as a useful diagnos¬ 
tic of non-thermal broadening, see Section 5.4. Note that the 
double-peaked line profile is caused by the velocity gradient 
and not by optically thick line formation. 

5. DIAGNOSTIC POTENTIAL 

In this section we present synthetic profiles calculated from 
the 3D atmospheric snapshot and how observables are related 
to physical quantities in the atmosphere, thereby exploring the 
diagnostic potential of the O I 135.56 nm line. Figure 7 sum¬ 
marizes the properties of the synthetic 01135.56 nm line cal¬ 
culated throughout the 3D simulation box by showing the to¬ 
tal intensity, the Doppler shift of the maximum emission and 
the total line-width (given as 1/e width, see Equation 11 for 
a definition). The total intensity is typically 2mW m~^ sr“^ 
with a typical range of 0.3 — 65 mW m~^ sr“^. The Doppler 
shift is up to 10 km s“^. The total line-width is over a large 
part of the simulation box 2 — 4 km s“^ with maximum width 
up to 12 km s“\ typically in the comers of the simulation box 
where the magnetic field is the weakest. These locations have 
acoustic shocks propagating through the line-forming region, 
causing large line-widths. 

Before we look at relations between observables and atmo¬ 
spheric properties, we inspect the contribution function to in¬ 
tensity to establish from where in the atmosphere observables 
are encoded. 

5.1. Contribution function to intensity 

Typically the O I line has a single emission peak without 
central reversal. We define the line core to be the location of 
maximum emission. Figure 8 shows the contribution func¬ 
tion to intensity (Equation 7) at this core wavelength for a cut 
through the 3D atmosphere at v = 12 Mm (same cut as used 
in Rathore & Carlsson 2015; Rathore et al. 2015). 

It is obvious that the intensity in the core of the 
01135.56 nm line mostly comes from heights above the T = 1 
height. The OI line thus shows optically thin formation. 

Correlations with atmospheric conditions are often made 
with the conditions at optical depth unity. For the 
01135.56nm line, with its optically thin formation, this is ob¬ 
viously a bad choice. Instead we define the formation height 
from the first moment of the contribution function: 

/ oo 

zCi^{z)dz. (9) 

We correlate with mean atmospheric parameters calculated 
in the same fashion as a contribution function weighted aver¬ 
age: 

Vfm = r V{z)c,,{z)dz (10) 

J —oo 

where V (z) represent the physical quantity along a column 
of the atmosphere, e.g. the electron density, the temperature, 
or the vertical velocity, C/^(z) is the contribution function to 
intensity and Vf^ is the contribution function weighted aver¬ 
age of the physical quantity V. 
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Figure 7. Total intensity (left), Doppler shift of the maximum emission (middle) and total line-width (right) of the 01 135.56nm line. 



Figure 8. Contribution function to the intensity of the O 1 line core along a 
cut through the atmospheric model at X = 12 Mm together with the z — 1 
height of the line core (red solid line) and temperature contour at 30 kK 
(green). The contribution function has been individually scaled for each col¬ 
umn for increased visibility. Note that the Z-axis has been stretched. 

5.2. Velocity 

Figure 9 shows that there is a tight correlation between the 
Doppler shift of the line core (defined to be the wavelength of 
maximum emission) and the contribution function weighted 
vertical velocity. This makes the O I 135.56nm line a good 
velocity diagnostic. The shift of the O I line is usually small 
(90% of the data points lie within ± 5 km s“^) and the av¬ 
erage shift is 1.5 km s“^. This average shift is caused by a 
global oscillation that is of this value at the time of the sim¬ 
ulation snapshot. The temporal average profile over a num¬ 
ber of snapshots shows a shift very close to zero. Using the 
O I 135.56nm line as a wavelength reference, as is done for 
IRIS, is thus warranted. 

5.3. Electron density 

As discussed in Section 4, the main channel of the OI emis¬ 
sion is through cascading from the continuum. The OI ioniza¬ 
tion degree is proportional to the hydrogen ionization degree, 
hence the population of the OI continuum will be proportional 
to the proton density. In the chromosphere, hydrogen is the 
dominant electron donor and the proton density can be well 
approximated by the electron density. The radiative recom¬ 
bination introduces another electron density dependency. We 



Figure 9. Probability Density Function (PDF) of the vertical velocity as a 
function of the Doppler shift of the line core. Each column is scaled to max¬ 
imum contrast to increase visibility. The inner green contour encompasses 
50% of all points and the outer contour 90%. The Pearson linear correlation 
coefficient is given in the upper left corner. The red line denotes the line y=x. 

therefore anticipate the 01135.56 nm line total intensity to be 
proportional to the square of the electron density. Figure 10 
shows the PDF of the logarithm of the contribution function 
weighted electron density as a function of the logarithm of 
the O I total line intensity. From the linear regression we get 
a slope 0.51 in the log^gA/g = /(log^Q/) plot, equivalent to 
I ^ This confirms that the total intensity indeed is pro¬ 
portional to the electron density squared. There is, however, 
a substantial width to the relation because of the temperature 
dependency of the recombination coefficient, see Figure 2. 

5.4. Non-thermal broadening of the OI line profile 

























8 


Lin & Carlsson 



log^o( I [W m' sr' ]) 



Figure 10. PDF of the logarithm of the contribution function weighted elec¬ 
tron density as a function of the logarithm of the OI total line intensity. Each 
column is scaled to maximum contrast to increase visibility. The green con¬ 
tours encompass 50% and 90% of all points. The Pearson linear correlation 
coefficient is given in the upper left corner. The red line denotes the line 
logioA^e = 0.511ogiof + 12.10, which is equivalent to / ~ • For rnore 

details see discussion in Sec. 5.3. 

In Section 4.3 we saw that the O I line-width can be sen¬ 
sitive to the velocity amplitudes in the line formation region. 
Here we examine this relation further. 

The total broadening, Wtotah thermal broadening Wth, and 
non-thermal broadening W^th, are defined as follows: 


f^total — 




41n(2) ’ 

Wa, = ^/2kT/m, 


Wnth = A/H^.Lr 


-Wth\ 


( 11 ) 

( 12 ) 

(13) 


where k, m, T, Wfwhm denote the Boltzmann constant, the 
atomic mass, the temperature and the full-width at half maxi¬ 
mum from our O I line profile, respectively. 

The measure of the velocity profile V (z) that is relevant to 
the line formation is defined as follows: 


V^ = V2 


J{5V{z))^C{z)dz 

JC(z)dz 


where 


(14) 

(15) 

(16) 

where Vw is V2 times the weighted velocity standard devi¬ 
ation, Vij and are the blue- and the red-side frequency that 
defines the FWHM and Ci^{z), V{z) and Vcore are the contri¬ 
bution function to intensity, the vertical velocity in the atmo¬ 
sphere and the Doppler shift of the defined line core. For a 
Gaussian velocity distribution, V^v would give a 1/e broaden- 


C(z)= rC/,(z)Jv, 

JVb 

5y(z)=y(z)-Vcore, 


Figure 11. PDF of (Equation 14) as a function of the total broadening, 
Wtotai (Equation 11). The green contours encompass 50% and 90% of all 
points. The Pearson linear correlation coefficient is given in the upper left 
corner. The diagonal red straight line denotes y = x, which is the expected 
relation for a Gaussian velocity distribution; the blue straight line denotes y 
= 0.68 X, which is the expected relation for a flat velocity distribution. 



Figure 12. PDF of (Equation 14) as a function of the non-thermal broad¬ 
ening, Wnth, of the line profile (Equation 13). Each column is scaled to maxi¬ 
mum contrast to increase visibility. The green contours encompass 50% and 
90% of all points. The Pearson linear correlation coefficient is given in the 
upper left corner. The diagonal red straight line denotes y = x, which is the 
expected relation for a Gaussian velocity distribution; the blue straight line 
denotes y = 0.68 x, which is the expected relation for a flat velocity distribu¬ 
tion. 

ing of the absorption profile of the same amount, which makes 
this definition appropriate for a direct comparison with W^th as 
defined above. 

Figure 11 shows the PDF of the measure Vw as a function of 
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the total line-width of the intensity profile of the 01135.56 nm 
line. The relation expected for a Gaussian velocity distribu¬ 
tion without thermal broadening is shown with a straight red 
line. The relation expected for a fiat velocity distribution (all 
the velocities are equally probable, e.g., a linear velocity gra¬ 
dient) is shown with a blue straight line. It is clear that when 
the total line-width is small (< 6 km s“^), the total line-width 
is larger than the velocity amplitudes indicate because the 
broadening is dominated by thermal broadening rather than 
non-thermal broadening. As the line-width grows bigger the 
broadening gets more and more dominated by non-thermal 
broadening. 

With the knowledge of the temperature structure in the at¬ 
mosphere, we may directly correlate the velocity field with 
the non-thermal broadening (using Equation 13). Figure 12 
shows this correlation together with the expectations from 
a Gaussian and a fiat velocity distribution. The correlation 
from our simulation shows a tendency of a fiat distribution 
when the non-thermal broadening is small (<4 km s“^), and 
it gradually approaches a Gaussian distribution for larger non- 
thermal broadening. Thanks to the optically thin formation, 
the 01133.56 nm line-width is a very good diagnostic of non- 
thermal broadening. 

6. COMPARISON WITH OBSERVATIONS 

We have explored correlations between atmospheric param¬ 
eters in a snapshot of a radiation-MHD simulation of the so¬ 
lar atmosphere and observables of the O I 135.56nm line. 
A relevant question is how well the simulated line-profile 
reproduces the real, observed, intensity profile. Figure 13 
shows the average synthetic O I 135.56nm line profile com¬ 
pared with the average profile from an IRIS observation. The 
synthetic profile has been convolved with the IRIS spectral 
point-spread-function (FWHM of 26 mA, De Pontieu et al. 
2014). The IRIS data were acquired on 2014 August 10 
starting at 03:12:39 UT and comprise a 400-step raster of a 
rather quiet area close to disk center. The exposure time was 
30 s. The standard IRIS wavelength calibration is based on 
the O I 135.56nm line such that the mean Doppler-shift is 
zero by definition. The ’’continuum” signal is dominated by 
contamination light from longer wavelengths. We have there¬ 
fore subtracted this constant level from the observed profile 
and added the continuum level from the simulations. The 
synthetic profile is more narrow than the observed profile. 
This deficiency of the simulation is also seen in profiles of 
the Mg II k line (Feenaarts et al. 2013a,b; Pereira et al. 2013) 
and the C II multiplet at 133.5 nm (Rathore & Carlsson 2015; 
Rathore et al. 2015). Given the good correlation between the 
O I 135.56 nm line non-thermal width and the velocity field 
in the atmosphere this is an indication that the simulation has 
too small macroscopic velocities. 

7. CONCLUSION 

In this work we have studied the basic formation mecha¬ 
nism of the 01135.56 nm intersystem line. The line shows an 
optically thin formation and the emission gets contributions 
from the whole chromosphere. The radiative rate is domi¬ 
nated by a recombination cascade and the total intensity is 
therefore proportional to the electron density and the popu¬ 
lation density of ionized oxygen. The ionization balance is 
dominated by charge transfer collisions with hydrogen such 
that the final total intensity is proportional to the electron den¬ 
sity squared. 



Figure 13. Comparison of the average synthetic spectrum from our simula¬ 
tion convolved with the IRIS PSF (red) with observations of the quiet sun at 
disk center by IRIS from 2014 August 10 at 03:12:39 UT (black). 

The Doppler shift of the OI 135.56nm line is very closely 
related to the velocity averaged over the formation region. 

Due to the optically thin nature of the O I 135.56nm line 
formation we get a line-width that is a direct measure of the 
thermal and non-thermal broadening. The thermal width is 
rather small (3.2km s“^ at a temperature of 10 kK, which is 
a temperature at the upper range of what we find) such that a 
profile with an 1/e width of more than 5 km s“^ is dominated 
by non-thermal broadening. The O I 135.56nm line thus pro¬ 
vides a very good measure of non-thermal velocities in the 
chromosphere and is an excellent complement to the optically 
thick IRIS chromospheric diagnostics such as the Mg II h & k 
lines and the CII lines around 133.5 nm. 

The research leading to these results has received fund¬ 
ing from the European Research Council under the European 
Union’s Seventh Framework Programme (FP7/2007-2013) / 
ERG grant agreement no 291058. This research was sup¬ 
ported by the Research Council of Norway through the grant 
’’Solar Atmospheric Modelling” and through grants of com¬ 
puting time from the Programme for Supercomputing. 
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